Photon-Photon Correlations as a Probe of Vacuum Induced Coherence Effects 



Sumanta Das and G. S. Agarwal 
Department of Physics, Oklahoma State University, Stillwater, Oklahoma H018, USA 

(Dated: February 2, 2008) 

We present new experimental implications of the effects of vacuum induced coherence on the 
photon -photon correlation in the vr-polarized fluorescence in j = 1/2 to j = 1/2 transition. These 
effects should be thus observable in measurements of photon statistics in for example Hg and Ba 
ion traps. 
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I. INTRODUCTION 



An early work fT] had predicted a very unusual effect 
of quantum interference in the problem of spontaneous 
emission. It was for example shown that in a degenerate 
V-system one could get population trapping and gener- 
ation of quantum coherences in the excited states. One 
of the key conditions for the occurrence was that the 
dipole matrix elements of the two transitions from the 
excited states were orthogonal. The later condition is 
difficult to meet though very large body of theoretical 
literature has been devoted to the subject of vacuum 
induced coherences [1, H, 0, B • It was also suggested 
how the above condition on dipole matrix elements 
can be bypassed if we consider anisotropic vacuum Q 
which for example would be the case while considering 
emission from excited atoms on nano particles Q. It is 
clearly important to find out easily realizable systems 
so that experimental results can be obtained. Kiffner. 
et. al.0 showed that one very important case would 
involve j = 1/2 to j = 1/2 transition. They calculated 
the spectrum of the emitted radiation and showed how 
vacuum induced coherences change the spectrum of the 
emitted radiation. There are many systems where it 
is easy to find j = 1/2 to j — 1/2 transition. In fact 
in an early experiment of Eichmann et. al.fioj such 
transitions in ^^^Hg^ were used to examine interferences 
in a system of two ions. A recent work uses the j = 
1/2 to j = 1/2 transition in ^^^Ba+ ions [ll|. While 
the results of Kiffner. et. al. on the spectrum are quite 
interesting the current experimental efforts [ll|, [l^, 
are more focussed on the study of photon-photon 
correlations. Thus an important question would be - do 
the vacuum induced coherences significantly affect the 
photon-photon correlations ? This is the question we 
answer in affirmative. 

It may be added that the photon-photon correlations 
have acquired new significance in the context of quantum 
information processing and quantum imaging as well as 
in interferences from independent atoms [l3 . [isl . [l^ . 
Thus it is important to have a consistent calculation of 
such correlations in situations where vacuum induced 
coherence (VIC) effects are important. 
The organization of this paper is as foUows-In Sec 
2 we introduce the model and present the working 



equations. In Sec 3 we calculate the photon-photon 
correlations both in presence and in absence of the 
vacuum induced interference effects. In Sec 4 we present 
numerical results to highlight the effects of vacuum 
induced coherences on photon-photon correlations. In 
Sec 5 we conclude with the outlook and future directions. 



II. MODEL 



ll/2-l/2> = ll> 



I A I I1/2,1/2> = I2> 



CO 



2y 




ll/2-l/2> = l3> 



I1/2,1/2> = I4> 



FIG. 1: Schematic diagram of a four-Level atom modelled by 
j = l/2toj = l/2 transition 



The Fig.l. shows the level scheme of a four-level 
atom modelled by j = 1/2 to j = 1/2 transition. This 
kind of level scheme is easily realizable and has already 
been studied, for example in ^^^Hg+ ^ and "^Ba+[li| 
ions. The ground level is 6s^S'i/2 and the excited level 
is 6p^Pi/2- Each of these levels are two fold degener- 
ate. The transitions |1) ^ |4) and |2) ^ |3) couple to 
cr"*" and a~ polarized light respectively. The transitions 
|1) ^ |3) and |2) ^ |4) couple to fight linearly polarized 
along the and their dipole moments are antiparallel. 
The spontaneous decays of the excited state to the two 
ground states are given by 27 and 27^ as shown in the 
figure. The electric dipole moment operator for this level 
scheme is defined as 



A, 



(»,J = l,-,4) 



(1) 
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The non- vanishing matrix elements of the electric dipole 
moment operator d can be found using the Wigner- 
Eckart theorem and are given by, 



d42 = 



dil — d^2 = J 



V6 
1 

73 



(2) 



with e_ = {x — iy)/^/2. In Eq. (2) T> denotes the reduced 
matrix element of the dipole moment operator d. The 
four-level system is driven by a tt polarized monochro- 
matic field of frequency oj , 



E(<) = £oe"'"'e^ -I- c.c. , 



(3) 



were c.c is the complex conjugate. With this particular 
choice of polarization, the driving field couples only to the 
two antiparallel dipole moments d^i and c?42. The total 
Hamiltonian for this atom-field system is then given by 

n^HA + ni , (4) 

where the unperturbed Hamiltonian for the atom is. 



Ha = h^uji\i){i\ , 



(5) 



an the interaction Hamiltonian for this system is given 
by 



Hi = -d.E(t) 

= ftJ7,(|l)(3|-|2)(4|)e— * 



H.c 



(6) 



where H.c is the Hermitian conjugate and f2c is the Rabi 
frequency defined by 



£o 



(7) 



The time evolution of this four level system is investi- 
gated by studying the density matrix equation. The 
spontaneous emission is included via the master equa- 
tion techniques. Following the standard procedure [l| we 
obtain, 



elements di3 and ^24 are anti-parallel. In a frame rotat- 
ing with the frequency of the coherent drive the density 
matrix equations are, 



Pii = i^*cPi3 - i^cPsi ~ 2r(5ii , 

PZZ = i^cPZl - i^*cPl3 + laP22 + 7Pll , 

Pl2 = ~?ricP32 - i^*cPlA " 2rpi2 , 

/5i3 = -zApi3 -f zric(pii - Pas) - rpi3 , 

pi4 = -iA/5i4 - i^ch2 - i^cPsi - r/5i4 

P23 = -i^P2Z + i^chl + i^chz - rp23 

P24 = ~i^P2A - i^c{P22 - Pii) - rp24 , 

P34 = ~i^cP32 - i^*cPl4 - 79P12 , 



(9) 



where 



Pii — Pii 1 Pij 

r = (70- + A^UJ- CJ13 : 



— iujt 
UJ — UJ24 , 



(10) 



The remaining equations can be generated by taking 
complex conjugates and using Tr{p} = 1. The steady 
state solution of Eq. (10) are found to be 



P12 
Pii 

P33 
Pl3 



Pii = P32 = P3i = , 

1 

P22 



(11) 



2[2|r!e|2 + r2 + A: 



Pii = n 



2[2|r!e|2 



-P2i = 



r2 



(12) 



r2 



r + iA'2[2|i7c|2 + r2 + A^ 



As can be seen from Eqs. (11) and (12) the vacuum 
induced interference has no effect on the steady state 
solutions. Clearly vacuum induced coherence can show 
up in dynamical quantities. 



III. PHOTON-PHOTON CORRELATIONS 



i 

P = -jXKiP]+i^P : 

Cp = -7,[|l)(l|p+|2)(2|p + p|l>(l| 

-f p|2)(2|-2|3)(3|p22-2|4)(4|pn] 

- 7[|l)(l|p+|2)(2|p + p|l)(l| 

+ p|2)(2|-2|3)(3|pn-2|4)(4|p22] 

+ 7[|4)(3|p2i + |3)(4|pi2] , (8) 

The last two terms in Eq. (8) arise from the vacuum 
induced interference and it comes as the dipole matrix 



Since the objective of this paper is to investigate the 
observable consequences of the vacuum induced coher- 
ence ; we focus our attention on the photon statistics of 
the radiation emitted by our model system. We in par- 
ticular will calculate photon-photon correlations as cur- 
rently considerable experimental effort is on such correla- 
tions. For this we need to know how to relate the atomic 
properties with the statistical properties of the sponta- 
neously emitted radiation. The answer to this question 
already exists in quantum theory. In fact from the exist- 
ing [l5| literature, we know that the positive frequency 
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part of the electric-field operator at a point r in the far- 
field zone can be written in terms of the atomic operators 
as 



E+(r,t) 



(r,t) - klY^{\R, X {R, X cTsOvlai] 



+ \Ri X {Ri X ^42)^42] 

+ {Ri X (Ri X ^32)^32] 

+ [R^ X {R, X d4l)A4i]}R-^ 

V p-iikaf-fi+uT) 



(13) 



where Ri — f ^ fi , f being the distance of the point of 
observation from the origin and fi being the position of 
the atom from the origin. Further t — t — ^ is the re- 
tarded time, ko = luq = W13 — (^24, dij is the electric 
dipole moment operator and the atomic operators are as 
defined in Eq. (1). The first term on the right of Eq. (fT3|) 
is the free field term and the second term is the retarded 
dipole field emitted by the atom. The emitted radiation 
consist of different polarization components- the n and 
the a polarized components. The terms A31 and ^42 
correspond to tt polarization whereas the ones ^32 and 
A41 correspond to a polarization. We next calculate the 
photon-photon correlations and the normalized second 
order correlations for the emitted radiations from the tt 
transitions of this driven four-level atom. For n polar- 
ization the relevant part of the electric field operator is 
given by, 



E+(f,i) = E+(f,<) - {"^r-iin X (n x d3i)]|3)(l| 



[nx (rixd42)]|4)(2|J 



(14) 



In the lowest order correlation the free field term of Eq. 
(fT4| does not contribute. This can be seen directly from 
the definition of quantized fields [l^, the fact that the 
field is initially in the vacuum state and the expression 
for the normally ordered correlation function for the field, 
(E^(r, t)E+(r', t')). Hence with no contribution from the 
free field term the intensity /^r of the light emitted on the 
TT transition from the atom is , 

(/.) = (E-(f,i).Ej(f,<)) 

= X (n X d3i)]* • [n X (n X d3i)]|l)(l|, 

c r 

+ [n X (n X ^42)]* -[fix (fix d42)]|2)(2|J , (15) 

where we have taken our origin at the location of the 
atom , r — fir, r is the retarded time and we used the 
property AijAki — AnSkj- The negative frequency part 
of the electric field operator E~ (r , t) can be found by 
taking the complex conjugate of the positive frequency 
part. Now if we assume that the point of observation lies 
perpendicular to both the polarization and propagation 
direction we have from Eq. (jlSp 



Eq.(16) can be further simplified using Eqs. ^ and (|12p . 
where in using Eq. (|12p we have assumed that observa- 
tion is been made at long time limit. The final expression 
for in the long time limit (steady state) is then. 



6r2 [2|f7c|2 + r2-H A2] ' 



(17) 



(I^) = (-)'i(|rf31P(|l>(l|)r + M42P(|2)(2|)0 , (16) 



Eq.(17) clearly show that intensity emitted on the tt tran- 
sitions is not altered by vacuum induced coherences and 
is simply proportional to the steady state population of 
the excited states. 

Let us now investigate what happens incase of two time 
photon-photon correlations on the tt transitions. The 
two-time photon-photon correlation for the level scheme 
in Fig.l can be written as 

{U{t + T)hit)) = (E-(r,t)E~(f,t-Hr) : 
Ej(r,i + T)Ej(f,i)) 

= {-r^{[nxihxd,,)]* 
c 

■[h X (n X dsi)]}'^ 

((|1)(3|-|2)(4|),(|1>(1| + |2)(2|),+. 
(|3)(1|-|4)(2|),), (18) 

The two-time correlation function which appears in Eq. 
(18) is to be obtained from the solution of the time- 
dependent density matrix equations (Eq.(9)) and the 
quantum regression theorem .1(|. A closer look at Eq. 
(9) show that eight of the fifteen equations form a 
closed set of linear equations which can be solved to find 
|l)(l|t+r, |2)(2|t+r and hence the term (|l)(l|-h|2)(2|)t+^ 
in Eq. (18). Before going further let us list those eight 
equations, 

Pii = i^*cPi3 - i^cPsi - 2rpii , 

P33 = i^cPai - i^lh'i + 7<tP22 + 7Pii , 

/5i3 = -zA/5i3 -I- zric(/5ii - P33) - r/5i3 , 

P31 = iAp3i - zri*(pii - P33) - r(53i , (19) 

P22 = i^cPi2 - i^*cP2i - 2rp22 , 

P44 = iil*P24 - i^cPi2 + laPll + 7P22 , 

P24 = -iAp24 - i^c{p22 - P44) - rp24 , 

P42 = iAp42 -I- ifl'*{p22 - P44) - rp42 , 

In compact notation this equations can be written as, 

P = Mp, (20) 

where p , p are (8x1) column matrix and A4 is a 
(8x8) square matrix. Now using the method depicted 
in [17] and using Eg. p^ the solution of (|l)(l|t-|_^) and 
{\2){2\t+T) can be expressed in the form 

= /ii(r)(|l)(lk)+/i2(r)(|3)(|3|*) 

+ /i3(r)(|3)(l|,) + /i4(r)(|l)(3|,) 

+ A5(r)(|2)(2|,)-H/i6(r)(|4)(4|,) 

+ /iv(r)(|4)(2|,)-H/i8(r)(|2)(4|,),(21) 
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(|2)(2|,+0 = /5i(r)(|l)(l|,) 
+ /53(r)(|3)(l|,) 
+ /55(r)(|2)(2|i) 
+ /57(r)(|4)(2|,) 
where the f 's are defined by 

Mr) - (e^- 



■/52(r)(|3)(|3|0 
■/54(r)(|l)(3|t) 
■/56(r)(|4)(4|,) 
■/58(r)(|2)(4|,),(22) 



where we have used Eq. (2) for the dipole matrix ele- 
ments. Note that F2{t)[Fq{t)] is the sum of probabilities 
of finding the atom in the states |1) and |2) given that at 
T = 0, the atom was in the state |3)[|4)]. In the limit of 
large r, 



and 



)ik : 



^Ik 



(23) 



(24) 



G(2)(r)-^(ii^)8&(- 



r2 



(29) 



Here we have diagonalized the matrix A4 with A being 
the eigenvalues and P being the corresponding eigenvec- 
tors. We now make use of the quantum regression theo- 
rem to obtain the two time correlation function as, 

{BHtm(l\ + \2)(2\UrB{t)) ^ F,{r){BHt)\ 

+ F2iT){BHt)\ 

+ F3iT){B^{t)\ 

+ F,{r){BHt)\ 

+ F,{T){BHt)\ 

+ FeiT){BHt)\ 

+ Fr{T){B\t)\ 

+ FsiT){BHt)\ 



Next let us derive the expression for two-time photon- 
photon correlation in absence of interference. In this case 
the total photon-photon correlation will be a simple addi- 
tion of photon-photon correlations for radiation emitted 
on individual tt transitions. 



|3)(3|ti3(t)) 
|3)(l|tB(i)) 
\l){i\tB(t)) 
|2)(2|tS(t)) 
|4)(4|tB(i)) 
|4)(2|tS(i)) 
|2)(4|*S(i)) , 
(25) 

where we define the operator B as, B^{t) = (|1)(3| — 
|2)(4|)t ; B{t) = (i?t(t))t and F,(r) = fu{r) + Mr). 
Using this new definition of the operator in Eq. (18), 
the expression for the two-time photon-photon correla- 
tion becomes, 

(/.(i + r)/.(t)) = {^)^^{[nx{hxds,)]* 
c 

■[n X {n X d3i)]f (26) 
X {B\mi){l\ + \2){2\),+rB{t)) , 

which when Eq. (25) is used, simplifies to 

(/,(< + r)/,(0) = (i^)8l{[r^ X (fi X dsi)]* 
c 

■[hx{nx d3i)]}2 (27) 
X (F2(r)(|l)(l|),+F6(r)(|2)(2|),), 

In the long time limit (|l)(l|)t = pii(i) and (|2)(2|)t = 
p22{t), where pii{t), p22{t) are the steady state popula- 
tions of the excited states given by Eq. (fT2|) . Now fol- 
lowing our assumption that the point of observation lies 
perpendicular to both the polarization and propagation 
directions and substituting for pn , p22 from Eq. (|12p . 
Eq. (27) can be further simplified. The final expression 
for the two-time photon-photon correlation is then, 

Gi'Hr) = {lAt + T)h{t)) 

- (-)^S(^2(r) + Fe(r)) 
c 3Dr* 



G(2)(r) ^ {I^{t + T)I^{t)) 

= (E-(f,t)E-(f,t + r) : 

E+(f,i + T)Ej(f,i))|i)<3| 
+ (E-(r,t)E-(r,t + r) : 

E+(f,i-fT)E+(f,i))|2>(4| 



= (^f\{[nx{nxd,,)r 
c 

■[n X {h X dsi)]}^ 
(|1)(3|,(|1)(1|),+.|3)(1|,) 

+ (^)«^{[nx(nxd42)]' 

■[hx {nx d42)]V 
(|2)(4|,(|2)(2|),+.|4)(2|,) 



(30) 



(31) 



Finally using Eq. (21), (22) and (12) we get the photon- 
photon correlation in absence of interference as 



G?'(r) = (^^)«f^(/i2(r) + /56(r)) 
c 6or^ 



X ( 



2[2|r!e|2 + r2 + A2] 



(32) 



Here /i2(t)[/56(t)] is the probability of finding the atom 
in the states |1) [|2)] given that at t = 0, the atom was 
in the state |3) [|4)]. Eq. (32) in the limit of large r 
becomes. 



K'[r} (— ) — — ( 



36r4^[2|f^,|2 + r2 + A2] 



(33) 



We now further calculate the normalized photon-photon 
correlation corresponding to Eq. (28) and Eq. (32). The 
(7^2) function gives the non-classical aspects of photon 
statistics. 



1 



2[2|J7c|2-hr2 + A2] 



(28) 
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{I,{t + T)h{t)) {F2{t) + F,{t))p 
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{I^it + T)){lAt)) 

r 



(34) 



,(2) 



{t + T,t) 



{I^{t + T)I^{t)\ 



{{I^it + T)){lAt)))\l)i3l + {{I. 



Here pn is the steady state population of the excited 
state given by Eq. (12) and g^^^ [g^^''] is the normahzed 
two time photon-photon correlation function correspond- 
ing to presence [absence] of vacuum induced interference. 



IV. NUMERICAL RESULTS 

In this section we present our numerical results and 
discuss their consequences. To begin with, we first dis- 
cuss our method of computation. The decay rates of the 
excited states to the two ground states, and 27 are 
proportional to |d4ip and |c?3ip respectively. From Eq. 
(2) we get, 27cr = 7o/3 and 27 = 70/6 , where 70 is 
proportional to the square of the reduced dipole matrix 
element. We use these values for the decays in our nu- 
merical computation and normalize all the computational 
parameters with respect to 70. Further we use standard 
subroutines to diagonalize the complex general matrix 
A4 and obtain complex eigenvalues and eigenvectors of 
the form {a + i[3). For all values of detuning and Rabi 
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25 



30 



FIG. 2: Plot of two-time Photon-Photon correlation as a func- 
tion of time for Q. = 0.5, A = 0.0. All the parameters are nor- 

malized with respect to 70, where 70 — The blue and 

red lines correspond to photon-photon correlations in pres- 
ence and absence of VIC respectively. 



TABLE I: Eigenvalues for the diagonalized matrix M corre- 
sponding to two different values of the Rabi frequency of the 
driving field which is on resonance with the atomic transi- 
tions. 



A 


!^ = 0.570 


!^ = 3.O70 


1 
2 
3 
4 
5 
6 
7 
8 


(-0.349797,-1.10904) 
(-0.349797,1.10904) 
(-0.215794,-1.09726) 
(-0.215794,1.09726) 
(-0.300406,0.000000) 
(-0.165314,0.000000) 
(-0.403098,0.000000) 
(0.000000,0.000000) 


(-0.375000,5.99870) 
(-0.375000,-5.99870) 
(-0.208269,5.99522) 
(-0.208269,-5.99522) 
(-0.250000,0.000000) 
(-0.250000,0.000000) 
(-0.333462,0.000000) 
(0.000000,0.000000) 



frequency used in our computation we have two pairs of 
complex conjugate eigenvalues and four other eigenval- 
ues whose complex part are so small compare to the real 
part that these complex parts have no significant contri- 
butions. Hence these four eigenvalues can be taken to be 
purely real. Note that this is in contrast to the case of 



photon-photon correlations for the two level model where 
the number of eigenvalues is four [18]. The changes in the 
eigenvalues lead to spectral modification as discussed by 
Kiffner. et. al. The eigenvalues for VL^ = O.570, and 
= 370 and detuning A = are hsted in the Table (I). 
Note for example that for ~ 37o we have eigenvalues 
±5.99870i - 0.375 and ±5.99522i - 0.208269. This dif- 
ference in the real parts can produce a dip in the side 
bands in the MoUow spectrum. Next we calculate the 
elements fij of the 8x8 matrix [/] using Eq. (23) and 
Eq. (24). Finally we use the elements fij corresponding 
to Eqs. (28), (32) and Eqs. (34), (35) to evaluate the two 
time photon-photon correlations and normalized photon- 
photon correlations in presence and absence of vacuum 
induced interference respectively. 

The Figs. (2-4) show photon-photon correlations corre- 
sponding to Eqs. (28) and (32). The blue and red line 
in the figures correspond respectively, to photon-photon 
correlations in presence and absence of interference. The 
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FIG. 5: Probability for finding the atom in state |1) (/12) and 
|2) (/52) given that at time r = the atom was in the state 
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FIG. 4: Plot of two-time Photon-Photon correlation as a func 
tion of time for SI = 3.0, A = 0.0. All the parameters an 



normalized with respect to 70, where 70 — 



FIG. 6: Normalized photon-photon correlations plotted as a 
function of time, for Q — 0.5, A — 0.0. All the parameters are 

normalized with respect to 70, where 70 — ^^-^-x-^. 



correlations calculated in presence of interference shov 
strong damping of the oscillations and attain an over 
all higher value as the time separation r between twc 
counts increases. The differences between G^^^ and G'^- 
are most noticeable in the limit of large time separation 
T. In order to understand this we examine the distinction 
between F2{t) = fuir) + /52(t) and fuir). We recall 
that /i2[/52] was the probability of finding the atom in 
the state |1)[|2)] given that at r = 0, it was in the state 
|3). We exhibit these probabilities in the Fig. (5). We 
observe that the function fh2{T) starts becoming signifi- 
cant at the time scale of the order of 7"^. 



Further for large r, /12 and /52 become comparable. The 
physical process that contributes to f^2 is the following, 

TT—pol emission n—pol 

Similarly population can start from the state |4) and end 
up in the state |1) via, 



1 i In spr , XT—photo7i , 

4 2 3 



TT—pol emission ir—pol 



ID 



We show normalized photon-photon correlations in a typ- 
ical case in the Fig. (6) . In case of interference we observe 
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stronger damping of the oscillations and an overall reduc- 
tion of the g'-^-' function at shorter time scales. At long 
time limits 5^^-'(t — > oo) is 1. Photon antibunching effect 
is also visible as ^ g^^H^) < 1- For shorter time scale 
we get 5^^-'(t) ^ 1 a clear signature of the nonclassical 
nature of the two-time photon-photon correlations. 

V. CONCLUSIONS 

In conclusion we have shown that the vacuum in- 
duced coherence (VIC) do significantly affect the two-time 
photon-photon correlations even though they show no ef- 
fect on the total steady state intensity of the radiation 
emitted on the tt transitions. The effect of this coherence 



is reflected in form of stronger damping and overall larger 
values of the correlation function G*-^-* . The level scheme 
j = — 1/2 is easily realizable and has already 

been used, for example in ^^^Hg~^ [lO^ in the context 
of interferences produced by a system of two ions and 
more recently in ^^^Ba~^ [11] in the context of emission 
in presence of a mirror. In future we hope to investigate 
how the asymmetry in the level structures introduced by 
a magnetic field [19| would influence the photon-photon 
correlations. This might in turn give us more freedom 
in choosing the level structure and hence more broader 
choice in selecting atomic transition for experiments. Fi- 
nally note that it would also be interesting to examine 
the VIC effects in the context of nonlinear optical exper- 
iments using j=l/2 to j=l/2 transitions. 
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